The magnetic form factor of iron in SrFe 2 As 2 
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The iron magnetic form factor in SrFe2As2 has been determined by neutron diffraction and by 
density functional theory (DFT). As noted previously, the magnitude of the calculated moment using 
DFT is sensitive to the Fe-As distance. However, the shape of the calculated form factor is practically 
insensitive to the Fe-As distance, and further we show that the form factor closely resembles that 
of bcc iron, and agrees well with experiment. The spin density exhibits some anisotropy due to 
geometry and As hybridization. 

PACS numbers: 75.25,+z, 74.70.Dd, 75.50.Ee, 71.20.-b 



The magnetism and nature of the iron-3d electronic 
states in FeAs based pnictide superconducting materi- 
als is of interest. The states near the Fermi level are 
involved in the superconductivity and are also responsi- 
ble for the sensitivity of the magnetic moment to small 
structural changes. At first glance the proximity of the 
superconducting states to an antiferromagnetic (AF) or- 
dered state suggests a similarity to cuprates; however 
the Mott-Hubbard AF insulator state in the undoped 
cuprates is not amenable to density functional theory 
(DFT) approaches 1 2 , whereas DFT has been remarkably 
helpful in elucidating the electronic and magnetic behav- 
ior of FeAs related materials. As one example, we note 
that DFT calculations for the ambient and high pressure 
phase of CaFe2As2 predicted the disappearance of the 
magnetic moment in the collapsed phase!™ Mazin et al. 
remark in their study of LaFeAsOi-^F^, "magnetism in 
this compound is very itinerant, making the calculated 
magnetic energies and moments extremely sensitive to 
the approximations used and to fine details of the crystal 
structure". They recommend the LDA functional and 
using the theoretical As positions for "best results" P In 
this work we investigate the spatial extent of spin density 
by both measuring the magnetic form factor and calcu- 
lating it using the experimental and also the theoretical 
As positions. 

The orthorhombic crystal structure and magnetic or- 
dering found in SrFe2As2 is illustrated Fig. [T] Above 
~200 K the crystal is tetragonal, but undergoes a simul- 
taneous structural and AF transition below this temper- 
ature. The experimental value of the magnetic moment 
on each iron site in the AF state of SrFe2As2 has been 
estimated using neutron powder diffraction to be 1.01(3) 
and 0.94(4) fi^. There is some discussion if the 
AF/structure transition is first or second order 89 but it 
is not relevant to our measurements taken at low tem- 
peratures (~10 K). 

A high-quality SrFe2As2 single crystal synthesized by 
the FeAs flux growth technique^ was selected for this 
study. This crystal was also used in a recent neutron 
diffraction study of the nature of the magnetic and struc- 
tural phase transitions in SrFe2As2p^ The crystallinity 
and purity were characterized by Laue back-scattering, 



X-ray powder diffraction, magnetization and resistivity 
measurements. The mosaic is 0.29(1)° full width at half 
maxima for the (0 8) reflection!^ The elastic neutron 
scattering measurements were carried out on the HB-1A 
fixed-incident-energy (14.6 meV) triple- axis spectrome- 
ter using a double pyrolytic graphite (PG) monochroma- 
tor (located at the High Flux Isotope Reactor, HFIR, at 
the Oak Ridge National Laboratory, USA). Two highly 
oriented PG filters, one after each monochromator, were 
used to reduce the A/2 contamination. The beam col- 
limation throughout the experiment was kept at 48'- 
48 / -sample-40 / -136 / . The single crystal (~ 15 mg) was 
wrapped in Al foil and sealed in a He-filled Al can which 
was then loaded on the cold tip of a closed cycle refriger- 




FIG. 1: (Color online) SrFe2As2 orthorhombic crystal struc- 
ture: the green spheres are Sr atoms, yellow are As atoms 
and both blue(Fel) and red(Fe2) are Fe atoms but represent 
different magnetic moment directions. The gray plane is the 
one that was chosen for the spin density contour plot of Fig. 
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ator with (H L) in the Fmmm symmetry as the scat- 
tering plane. 

The magnetic form factor F(K) is a Fourier transfor- 
mation of spin density at a given K vector and is given 
as 
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where K is a reciprocal lattice vector which is equivalent 
to a momentum transfer vector of scattering experiments. 

We used the spin polarized linear augmented plane 
wave (LAPW) method 12 to obtain the electronic struc- 
ture and spin density. The Fourier transformation of 
the plane wave expanded spin density in the intersti- 
tial region was added to the form factor for the angular 
decomposed spin density inside the muffin-tin spheres. 
To calculate the self-consistent spin density, we used 
a full-potential LAPW method 13 with a local density 
functional.— The quantity RMT-Kmax, where Rmt is the 
smallest muffin-tin radius and K max is the plane-wave 
cutoff, was set to 8.0 and 500 k-points in the irreducible 
Brillouin zone were generated for integration using the 
tetrahedron method.^ We used muffin-tin radii of 2.3, 
2.1 and 2.1 a.u. for the respective Sr, Fe and As spheres 
and 0.01 mRy/cell of the total energy as the convergence 
criterion. The structural parameters were taken from re- 
ported X-ray powder diffraction results. 16 Since it is well 
known that calculated ma gnetic properties are very sen- 
sitive to the As positon^ 17 * 18 1 , we performed the magnetic 
form factor calculation with both the experimental and 
the energy minimized "optimized" As position. To ob- 
tain the optimized As position, As atoms were relaxed 
along ^-direction until the total force on As atom was 
less than 0.1 mRy/o.w. . 

To obtain the form factor of iron in SrFe2As2 we col- 
lected the integrated intensities of the rocking curves of 
nuclear and magnetic Bragg reflections to the highest or- 
der possible with the HB-1A spectrometer. The nuclear 
reflections, necessary to obtain the magnetic moment, 
were chosen to be as close as possible in the scattering 
angle 26 to that of the corresponding magnetic Bragg 
reflections. This reduces any errors due to geometry cor- 
rections, Debye- Waller factor (DWF), absorption effects, 
etc. The integrated intensity of a rocking curve of a nu- 
clear Bragg reflection at a reciprocal lattice vector q in 
a crystal is given by: 
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where V is the scattering volume of the crystal, is the 
unit-cell volume, 3>o(0) is the beam flux at the angle #, 
-Fn(q) is the structure factor, A is the beam wavelength, \i 
is the absorption length, sin(2#) is the Lorentz factor for 
a rotating crystal, and 2W = q 2 (?iq) 2 is the DWF. We 
collect all the constants and unknown q-dependent fac- 
tors in C(q), including the DWF. The integrated inten- 
sities of the rocking curves of the chosen nuclear Bragg 



reflections were normalized to the corresponding struc- 
ture factors and the Lorentz factors. The data were also 
corrected for the fact that the rocking curves of nuclear 
Bragg peaks in our measurements include the twinned (ft 
0) and (0 k 0) orthorhombic domains, for normalization 
of the magnetic reflections which are due to the (ft 0) 
domain only. We assume C(q) is a DWF-like function 
(namely, a gaussian) and fit the data by the nonlinear 
least square technique to obtain a smooth C(q) func- 
tion. The refined parameter (uq) in this procedure is in 
good agreement with the corresponding values obtained 
in a powder diffraction study of SrFe2As2P 



TABLE I: The q values, integrated intensities, and the corre- 
sponding magnetic structure factors at 10 K of the magnetic 
reflections observed in single-crystal SrFe2As2 with the AF 
structure as shown Fig. [I] 



Reflection 


q (A" 1 ) 


Intensity 


|F M | 2 sin 2 a; 


Form Factor 


(101) 


1.244 


128.6±3 


10.92 


0.9170±0.0096 


(103) 


1.9084 


242. 1±4 


41.53 


0.8576±0.0170 


(105) 


2.798 


110.9±3 


53.6 


0.7085=b0.0328 


(107) 


3.756 


30.5±1 


58.2 


0.4954±0.0632 


(109) 


4.7417 


7.91=b0.58 


60.38 


0.3204±0.1058 


(301) 


3.441 


1.19 ±0.38 


1.4 


0.5678±0.0512 


(303) 


3.733 


6.57±0.39 


10.9 


0.5273±0.0623 


(305) 


4.257 


5.28=b0.28 


23.3 


0.3799=b0.0860 



Similarly, the scattered intensity of a magnetic Bragg 
reflection can be expressed as: 
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where f(q) is the magnetic form factor at the mag- 
netic reciprocal lattice q, -F M (q) is the magnetic struc- 
ture factor^, and sin 2 a = 1 — (q • ft) 2 where q and ft are 
the unit vectors along the scattering vector (reciprocal 
vector) and the direction of the moment, respectively (in 
this case ft = (1 0)). In Table [l| we list all the ob- 
served magnetic reflections, their integrated intensities, 
and the corresponding magnetic structure factors taking 
the magnetic structure as shown in Fig. [I] Using these 
parameters in Table [I] and Eq. Q we calculated the value 
of the form factor gSf(\q\) as shown below in Fig. |4^b) 
(blue triangles). 

To estimate the magnetic moment we fit the measured 
form factor with Fe 2+ ionic form factor to extrapolate 
the data to q=0 yielding M= 1.07(2) /i B . We also used 
the shape of the calculated form factor (described below) 
to extrapolate the data yielding a magnetic moment of 
M= 1.04(1) fi B . 

We now turn to the theoretical calculations of the spin 
density and magnetic form factor. The theoretical energy 
minimization process of relaxing the As position resulted 
in zas = 0.3499 while the starting experimental position 
was zas = 0.3612. The distance between Fe-As was de- 
creased about 3% (from 2.392 Ato 2.315 A) and the mag- 
netic moment on the Fe atoms was changed from 1.676 
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/jLb to 0.879 /x B (at the optimized theoretical minimum 
energy za s position). Table |TT| shows the Fe ^-orbital 
decomposed charges and moments. Although there are 
variations among the occupied electronic orbit als, mag- 
netic moments are fairly well distributed across all of the 
orbitals. This would not be expected if one or more of 
the orbitals were "localized" . Indeed, all of orbitals con- 
tribute to a broad density of states. 

TABLE II: Fe d-orbital electron occupation and magnetic mo- 
ment orbital decomposition. 





spin up 
(electrons) 


spin dn 
(electrons) 


Magnetic Moment 
Ob) 




exp za s 


Opt ZAs 


exp za s 


opt ZAs 


exp za s 


Opt Z As 


d 


3.810 


3.450 


2.142 


2.574 


1.668 


0.876 


d z 2 


0.719 


0.618 


0.495 


0.513 


0.224 


0.105 


d x 2- y 2 


0.757 


0.684 


0.402 


0.491 


0.355 


0.193 


d X y 


0.769 


0.700 


0.502 


0.568 


0.267 


0.132 


d xz 


0.767 


0.718 


0.354 


0.476 


0.413 


0.242 


dyz 


0.798 


0.730 


0.389 


0.526 


0.409 


0.204 



Figure [2] shows the spin density along Fe to As. It 
indicates that although the spin density is strongly con- 
centrated inside Fe atoms, the spin density inside the As 
muffin-tin sphere is affected by the surrounding Fe atoms 
through hybridization. With antiferromagnetic ordering 
(see Fig. [I]), the As atoms do not have a net moment. 
There are, however, spin distributions surrounding the 
As site that depend on the surrounding Fe atoms. The 
inset of Fig. [2] shows the radially weighted spin density 
inside the Fe muffin-tin sphere with a high density along 
the Fe-As direction and a significant drop in spin density 
along the z-axis. It clearly shows anisotropy of the spin 
density inside of the Fe muffin-tin sphere. 

Figure [3] shows spin density contours in the gray plane 
shown in Fig. [I] Since we are interested in the hybridiza- 
tion between Fe and As, we have chosen a plane that 
is determined by (1,1,0) and (0,0,1) vectors. It shows 
the antiferromagnetically ordered Fe atoms and the effect 
of the ordered Fe atoms on As atoms. The hybridiza- 
tion between Fe d electrons and As p electrons produces 
an anisotropic spin density distribution around the As 
atoms. 

Figure [4] a) shows the form factors that were calculated 
with the experimental As position (red filled squares) and 
the optimized As position (blue open circles). Data were 
normalized with the form factor at |if|=0 correspond- 
ing to the magnetic moment of the Fe atom. The ar- 
row points out the two scattering vectors with similar 
magnitude but different intensities - demonstrating an 
anisotropy increase for the "optimized" As position with 
a smaller magnetic moment. The K vector which has 
higher intensity is (5,0,1) and the lower intensity one is 
(3,0,9). Figure [I] b) includes a comparison between the 
experimental form factor of SrFe2As2, that of bcc irorP^, 
and the calculated one using the experimental As posi- 
tion. Even though SrFe2As2 has a much more compli- 
cated electronic structure than bcc Fe, their normalized 
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FIG. 2: (Color online) Spin density along the Fe-As direction 
as calculated from LDA with the experimental atomic posi- 
tion. Although As does not have a net magnetic moment in 
the antiferromagnetic structure, its spin density distribution 
depends on the magnetic ordering of Fe. The inset shows 
spin densities inside the Fe muffin-tin sphere along specific 
directions demonstrating substantial anisotropy. 
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FIG. 3: (Color online) 2-dimensional spin density plot indi- 
cating hybridization between Fe and As atoms and anisotropy. 
Since the spin density is strongly concentrated in the small 
radius inside Fe atoms(See Fig. [2]) a small spin density con- 
tour range was chosen to reveal the low spin density structure. 
The black line between Fe and As is the axis for Fig. [2] 
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FIG. 4: (Color online) a). Comparison of the form factors 
calculated with the experimental atomic positions (za s = 
0.3612) and theoretical energy optimized one (za s = 0.3499). 
The arrow points to scattering vectors that provide evidence 
for fairly strong anisotropy in the form factor, b). Experimen- 
tal form factors for SrFe2As2 and BCC Fe, and the calculated 
form factors with the experimental As position. 



form factors are quite similar, with differences showing up 



at larger and more difficult to measure scattering vectors. 
The similarity with bcc Fe is particularly evident in the 
total 3d charge density within the muffin-tin sphere. This 
is 5.95 electrons for the experimental atomic positions in 
SrFe2As2, 6.02 electrons for the corresponding calculated 
energy optimized As position and 5.83 electrons for bcc 
Fe calculated with the same muffin-tin radius. 

We have determined the magnetic form factor of Fe 
in SrFe2As2 by neutron diffraction experiments and de- 
duced the Fe magnetic moment as 1.04(l)/i B - We also 
calculated the magnetic form factor by first principles 
electronic structure methods using both the experimen- 
tal and optimized As positions. While the magnitude 
of the calculated magnetic moments strongly depends on 
the As position, the normalized magnetic form factors 
were remarkably similar with each other and also agreed 
well with the normalized magnetic form factors from the 
experiment. The largest difference between the two the- 
oretical calculations was the noticeable increase in the 
spin anisotropy for the optimized As position. The ex- 
periments were limited by spectrometer geometry and 
did not measure far enough in reciprocal space to access 
these contributions. 
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